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Accurate prediction of where and when typhoons (or named hurricanes which form over the North Atlantic 
Ocean) will make landfall is critical to protecting human lives and properties. Although the traditional 
method of typhoon track prediction based on the steering flow theory has been proven to be an effective way 
in most situations, it slipped up in some cases. Our analysis of the long-term Chinese typhoon records 
reveals that typhoons, especially super typhoons (those with maximum sustained surface winds of greater 
than 51 ms" 1 ), have a trend to make landfalls toward warmer land in China over the past 50 years (1960- 
2009). Numerical sensitivity experiments using an advanced atmospheric model further confirm this 
finding. Our finding suggests an alternative approach to predict the landfall tracks of the most devastating 
typhoons in the southeastern China. 

China is one of the countries hit most frequently and severely by typhoons, strong tropical cyclones with 
maximum sustained surface winds of greater than 33 ms -1 that form over the Northwest Pacific Ocean. 
On average, China sees 7 to 8 landfalls of typhoons per year, impacting a population of 250 millions with 
property values over 10 trillion dollars 1 . With the rapid increase of economy and population in the southeastern 
coast of China, the loss of property and life in these regions caused by the landfalls of typhoons may increase in the 
future. During the past few decades, observational networks and numerical models have been used to predict 
typhoon landfall tracks 2,3 , but large uncertainties remain. One of the major uncertainties is the prediction of where 
the typhoon will make landfall, which is essential for local decision makers to choose among proper emergency 
measures, such as evacuation of coastal residents to avoid loss of life. In this paper we report a new indicator to the 
landfall locations of typhoons, in particular super typhoons that pose the most severe threat to human lives and 
properties. 

Traditionally, typhoon tracks are predicted based on the steering flow theory, in which a typhoon is considered 
as an imbedded vortex in the height- averaged large-scale winds 4 " 10 . This scheme proves instrumental, producing 
reasonable forecasts of typhoon tracks in many cases. However, it has failed in several high profile cases 11 . In these 
cases, the landfalls of the super typhoons did not follow the steering flow, contradicting traditional wisdom. Other 
factors that control typhoon landfalls remain elusive. 

Typhoons, as well other tropical cyclones, are conceived over warm ocean surfaces where plenty of energy and 
water vapor are available 12 . Typhoons are intensified when moving over warmer sea surface, while weakened over 
colder sea surface 13 " 15 . Compared to the widely recognized ocean surface effects, little is known of what role the 
land surface plays in regulating typhoon tracks. Specifically, it is unclear if warmer land surface can influence the 
tracks of typhoons and hence their landfalls. In this paper, we report our analysis of the relationship between land 
surface temperature (LST) and typhoon landfalls in China during the past 50 years (1960-2009). The long-term 
records of land surface temperature and typhoon tracks are retrieved from the quality assured archive collected 
and maintained by the Chinese Meteorology Administration (CMA). 

Results 

Fig. 1 shows the tracks of super typhoons and LST two days before landfalls in southeastern China for the decades 
of 1960, 1970s, 1980s, 1990s, and 2000s. It is interesting to see that, for 1960s, 1980s and 2000s, the majority of 
super typhoons clustered their tracks and make landfalls at a warmer area with LST higher than 31°C, while the 
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Figure 1 | Tracks of super typhoons and composite land surface temperature (color shaded, unit: Celsius degree) two days before landfalls in the 
southeastern China for the 5 decades since 1960. (a) 1960s, (b) 1970s, (c) 1980s, (d) 1990s, (e) 2000s. 



tracks of super typhoons show a diverse pattern for 1970s and 1990s 
when LST is lower than 31°C. To illustrate more clearly the features 
of such a "warmward landfall", we examine the actual landfall tracks 
and the tracks predicted toward the highest land temperature points 
for all super typhoons during 1960-2009 (see Supplementary Fig. S2 
online) and show four of them in Fig. 2 which have a recurvature 
feature before landfalls (see more cases in Supplementary Fig. S2 
online). To make a comparison with the traditional track prediction 
method using steering flow theory, the direction of the correspond- 
ing steering flow is also shown in each panel. We define oci as the 
track deviation angle for the track predicted toward the highest land 
temperature point and a 2 as the one for the track predicted by the 
steering flow (see Supplementary Fig. SI online). The point of the 
highest LST is determined within a coastal area of 30°X30° sur- 
rounding the cyclone center. For these super typhoons, the values 
of oci are comparable to or even smaller than those of a 2 . The good 



agreement between the LST -predicted and actual landfall tracks sug- 
gests that, in addition to the environmental steering flow, the warmer 
land surface may be another important factor regulating the landfall 
tracks of tropical cyclones. 

How significant is the land temperature indicator for typhoon 
landfalls? We further compare the LST-predicted tracks to the 
observed tracks for super typhoons as well as the weaker cyclones 
during 1960-2009. Those typhoons that only landed on Taiwan or 
Hainan Islands are excluded in this analysis, due to the lack of LST 
data for these islands. We bin all typhoons into three categories based 

on wind speedj vj: typhoons (| vj: 32.7—41.4 m/s), severe typhoons 

(|vj:41.5~51 m/s), and super typhoons > 51 m/s). For each 

bin, these typhoons are further grouped based on their deviation 
angles, in which smaller deviation angle indicates better agreement 
between the LST-predicted tracks and the observed ones (see 
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Figure 2 | The land surface temperature (color shaded, unit: Celsius degree), 500-hPa geopotential height (dashed contour, unit: gpm), 
the observed tracks, and the deviation angles of predicted tracks by the land surface temperature and the steering flow at the time of two days before 
landfalls of the super typhoons, (a) Opal (1962). (b) Billie (1976). (c) Omar (1992). (d) Herb (1996). A 9-point smoothing was applied to the land surface 
temperature for 3 times. The blue solid square denotes the location of highest land surface temperature within an area of 30° X 30° surrounding the 
typhoon center at the time of two days before landfall, and the red empty dot denotes the location of landfall. The red arrow represents the steering 
flow. Oil an d oc 2 are me deviation angles of the predicted tracks by the land surface temperature and the steering flow, respectively (see Supplementary 
Fig. SI online). 



Supplementary Table SI online). Considering that the distance (L) 
between the typhoon position of two days before landing and the 
landing location may affect the calculated deviation angle, we further 
bin all typhoons into two categories based on the distance: L 



< 800 km and L > 800 km. The statistical results are presented in 
Table 1. For L < 800 km, the percentages of all typhoons with devi- 
ation angle a x < 15° and 15° <a x < 30° account for 33% and 22%, 
respectively; for L > 800 km, the corresponding percentages are 38% 



Table 1 Statistics on the accuracy of the land-temperature-predicted (steering-flow-predicted) landfalls compared against the observed 
tracks for different categories of typhoons making landfalls on the mainland of China during 1 960-2009. ai (or 0C2) is the deviation angle of 
the land-temperature-predicted (or the steering-flow-predicted) tracks compared against the observed tracks, and L the distance between 
the typhoon position at the time of two days before landing and the landing location 

Category Of typhoons (No. of events) ai (a 2 ) < 1 5° 1 5° < ai (a 2 ) ^ 30° 30° < ai (a 2 ) ^ 45° 45° < ai (a 2 ) ^ 90° 


L<800 km Typhoons (51) 


31% (33%) 


18% (28%) 


24% (18%) 


24% (20%) 


Severe typhoons (23) 


26% (17%) 


22% (35%) 


26% (26%) 


22% (17%) 


Super typhoon (30) 


40% (47%) 


27% (30%) 


20% (13%) 


1 3% (7%) 


Total (104) 


33% (33%) 


22% (31%) 


23% (19%) 


20% (15%) 


L > 800 km Typhoons (23) 


44% (35%) 


26% (39%) 


13% (17%) 


1 7% (9%) 


Severe typhoons (20) 


25% (30%) 


40% (55%) 


20% (15%) 


1 5% (0%) 


Super typhoons (29) 


45% (59%) 


31% (21%) 


7% (14%) 


1 7% (7%) 


Total (72) 


38% (41%) 


32% (38%) 


13% (15%) 


1 7% (5%) 
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Table 2 | The mean deviation angles a x (a 2 ) (unit: degree) of land- 
temperature-predicted (steering-flow-predicted) tracks compared 
against the observed tracks for different categories of typhoons 
making landfalls on the mainland of China during 1960-2009. 
These mean deviation angles a x (a 2 ) are averaged overall deviation 
angles ai (or 0C2) from 4 days to 1 day before the landfalls of all 
typhoons located at a given range of latitudes 

Category Of 

typhoons 10°-15° 15°-20° 20°-25° 25°-30° 



Typhoons 14.5(18.7) 29.8(26.2) 35.0(34.6) 23.7(28.1) 

Severe typhoons 14.7(16.9) 25.7(26.7) 36.5(25.3) 24.6(28.1) 

Super typhoons 11.2(14.9) 18.4(19.6) 23.0(18.5) 39.5(29.6) 

Total 13.8(17.3) 25.3(24.3) 30.8(26.1) 30.0(28.7) 



and 32%, respectively. For both L < 800 km and L > 800 km, super 
typhoons have the highest percentage for the smallest deviation angle 
(a x ^15°), indicating that super typhoons may move toward warmer 
land surface most likely. For comparison, the statistics for the devi- 
ation angle a 2 of steering-flow-predicted tracks are also given in 
Table 1 (indicated in parentheses). It is surprising to find that the 
accuracy of the LST -predicted tracks is comparable to that of the 
steering-flow-predicted tracks (also see Supplementary Table S2 
online). Table 2 shows the mean deviation angles of the predicted 
tracks by both the LST and the steering flow compared against the 
observed tracks for different categories of typhoons making landfalls 
at low- middle latitudes (10°N-30°N). These mean deviation angles 
are obtained by making an average over all deviation angles (ai or a 2 ) 
from 4 days to 1 day before the landfalls of all typhoons located at a 
given range of latitudes (with an interval of 5°) during 1960-2009. It 
is found that the mean value (ai ) of a x is comparable to that (a 2 ) of oc 2 
at each range of latitudes, and both a x and a 2 are smaller at low 
latitudes than in middle latitudes. For both the steering flow 
predictor and the LST predictor, the accuracy of track prediction 
for super typhoons appears to be higher than that for weaker ones 
except for the latitude range of 25° —30°. By examining the highest 
LST for super typhoons with a Y < 15° (with a total number of 25) 
during 1960-2009, it is found that there are 19 out of 25 with the LST 
higher than 3 1 °C, accounting for 76% of the total (see Supplementary 
Fig. S3 online). Therefore, based on Fig. 1 and Fig. S3 (as well as 
Fig. 5a), we may recommend the land surface temperature exceeding 
31°C as a track predictor for super typhoons. 

To investigate whether the LST predictor and the steering flow 
predictor are identical to each other, we calculate the percentages of 



the right (positive) or left (negative) deviation of predicted tracks by 
the two predictors for all typhoons during 1960-2009. The percen- 
tages are averaged from 4 days to 1 day before the landfall of each 
typhoon located at a given range of latitudes. The results are shown in 
Fig. 3. The LST-predicted tracks have a larger likelihood of deviating 
right at low latitudes than at middle latitudes while they tend to 
deviate left more likely at middle latitudes than at low latitudes. 
The situation for the tracks predicted by the steering flow is opposite 
to that for the tracks predicted by LST, i.e., they have a larger like- 
lihood of deviating right at middle latitudes than at low latitudes 
while they tend to deviate left more likely at low latitudes than at 
middle latitudes. The reason may be primarily attributed to the 
meridional variation of the coriolis force, the so-called beta-effect, 
which is considered to have significant impacts on the movement of 
typhoons 16 " 23 . In the northern Hemisphere, barotropic modeling 
studies have shown that the beta- effect may make the cyclonic vortex 
move northwestward at a speed of 1-3 ms" 1 due to the asymmetric 
coriolis force exerting on the cyclonic circulation and the inward 
convergent flows 18 " 23 . Because the relative position of a typhoon to 
the land surface is dependent on the latitude where it locates, the 
predicted moving direction of a typhoon at low latitudes using the 
LST predictor may be north-northwest most likely, while that at 
middle latitudes may be west- northwest most likely. Therefore, the 
beta-effect may make the actual tracks toward the west (left) at low 
latitudes and toward the north (right) at middle latitudes, i.e., the 
beta-effect may make the LST-predicted tracks deviate to the right at 
lower latitudes and to the left at higher latitudes more likely com- 
pared to their actual tracks. For the steering-flow predictor, the dir- 
ection of the steering flow is also greatly dependent on the latitude: at 
low latitudes, the steering flow may mostly follow the westward flows 
at the southern edge of the subtropical high, thus the beta- effect may 
cause the actual tracks toward the north (right), i.e., make the steer- 
ing-flow-predicted tracks deviate to the left more likely compared to 
their actual ones; at middle latitudes, the steering flow may mostly 
follow the northward flows at the western edge of the subtropical 
high, thus the beta-effect may cause the actual tracks toward the west 
(left), i.e., make the steering-flow-predicted tracks deviate to the right 
more likely compared to their actual ones. These results imply that, 
although the two predictors of the steering flow and the land surface 
temperature are not totally independent to each other, they are not 
identical to each other either, which makes the land surface temper- 
ature be a valuable predictor for typhoon tracks supplementary to the 
steering flow. The beta-effect is important and should be considered 
deliberately no matter using the LST predictor or the steering flow 
one. 
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Figure 3 | The percentages of the right (positive) or left (negative) deviation of predicted tracks by the land surface temperature and the steering flow 
compared against the observed tracks for all typhoons during 1960-1009. The percentages are averaged from 4 days to 1 day before the landfall of each 
typhoon located at different ranges of latitudes (with an interval of 5°). 
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Figure 4 | The simulated tracks of Longwang (2005) from the control run (CNTR) and different sensitivity experiments Rl, R2, R3, LI, L2, and L3. (a) 

Initializing at 0000 UTC 30 Sept. 2005 (3 days before landing), (b) Initializing at 0000 UTC 1 Oct. 2005 (2 days before landing). Here Rl, R2, and R3 
denote the experiments in which the land surface temperature is increased by 5 K, 10 K and 15 K at the initial time on the right hand side of the observed 
landing location, respectively, while LI , L2 and L3 denote the experiments in which the land surface temperature is increased by 5 K, 10 K and 1 5 K on the 
left hand side, respectively (see Supplementary Fig. S4 online for more details). 



To further confirm our finding of the "warmward landing" trend 
of super typhoons, we conduct a set of numerical sensitivity experi- 
ments using an advanced atmospheric model called Weather and 
Research Forecasting (WRF) which includes comprehensive physical 
processes 24 . We choose the super typhoon Longwang (2005) for our 
sensitivity experiments. We increase LST by 5 K, 10 K and 15 K on 
the right hand side (denoted as Rl, R2 and R3) and the left hand side 
(denoted as LI, L2 and L3) of the observed landing location, respect- 
ively, at the time of 3 days and 2 days before the landfall of 
Landwang (2005), while keeping all the other circumstances the same 
as those of the control run (denoted as CNTR) (see Supplementary 
section 1.4 and Fig. S5 online). It is found that when increasing LST 
on the right (left), the typhoon track deviates to the right (left), and 
the more the LST is increased, the more the track deviates from 
CNTR (Fig. 4). The results of these numerical experiments confirm 
that super typhoons do have a trend to move toward the warmer land 
surface. 

What is the mechanism of the "wardwarm landing" for typhoons? 
It is well known that tropical cyclones tend to move toward warmer 
sea surface to feed on moisture and energy 1215 . We propose here that 
warmer land surface, with the intense moisture supplied by summer- 
time Asian Monsoon 25 , plays a similar role in attracting typhoons to 
move toward warmer ground surface. Fig. 5a shows the variation of 
the averaged land surface temperature 10 days before and after the 
landfalls of typhoons and super typhoons with a 1 <15°during 1960- 
2009. For super typhoons with ai<15°, it is found that LST starts to 
increase (up to 31°C) as early as 4 days before landfall, but drops 
sharply near landfall and increases again after landfall. The drop of 
LST could be due to the latent heat release by the strong convection 
and the cooling by rainfall on the ground during the typhoon land- 
fall. For typhoons, similar variation of LST is found but the increase 



of LST before landfall is relatively weaker than that for super 
typhoons, which indicates that the signal of LST as a track predictor 
for a typhoon is weaker than that for a super typhoon. One reason for 
this could be that the super typhoons, which contain tremendous 
moisture and energy, may have more intensified and farther-reach- 
ing cyclonic flows than weaker ones. These cyclonic flows could carry 
large amount of moisture to the warm land surface and thus favor the 
development of deep convection over the warm land (indicated by 
the large yellow arrows in Fig. 6). Therefore, the thermodynamical 
interactions between super typhoons and the warmer land surface 
could be more active and efficient, which makes the super typhoons 
move toward the warmer land surface more likely. Another reason 
could be that the vertical structure of the atmosphere over land is 
more unstable during super typhoons than that during typhoons, 
resulting in more deep convection over the warmer land surface. 
To verify this, we analyzed the temporal variation of the mean air 
temperature and water vapor anomalies during the landfalls of the 
super typhoons and typhoons with a!<15° (see Supplementary Fig. 
S7 online) as well as their differences (Fig. 5b) averaged over an area 
of 600 km X 600 km centered at the warmest land surface two days 
before landfalls during 2000-2009. Relatively wetter and colder air is 
found over the land surface several days before the landfalls of super 
typhoons with ai< 15°, compared to those of typhoons with ai< 15°. 
Since the mean land surface temperature before the landfalls of super 
typhoons is generally higher than that of typhoons, larger soil-air 
temperature difference occurred with colder air over the warmer 
land surface (see Supplementary Fig. S19 online), which, accompan- 
ied with more moisture probably supplied by both the southwesterly 
East Asian Summer Monsoon (EASM) and the detrained flows at the 
outer edge of super typhoons, may increase the instability of the 
atmospheric planetary boundary layer (PBL) of the atmosphere. 
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Figure 5 | Temporal variations of mean land surface temperature, air temperature and moisture averaged over an area with radius of 600 km centered 
at the warmest land surface two days before landfall, (a) Mean land surface temperature (unit: Celsius degree) 10 days before and after the landfalls of 
typhoons (black curve) and super typhoons (red curve) with a!<15° during 1960-2009. (b) Time-height cross section of the differences of the air 
temperature (color shaded, unit: K) and moisture (contour, unit: g'kg -1 ) between the super typhoons and typhoons with a!<15° during 2000-2009. 



Such an unstable structure of the atmospheric PBL favors the occur- 
rence of deep convection, which in turn feed the super typhoons with 
energy and moisture and thus "attract" the movement of the super 
typhoons. A diagnose of divergence and vorticity fields for the super 
typhoon Longwang (2005) using reanalysis data shows that the con- 
vergence center with positive vorticity of the super typhoon moved 
toward the convergence zone with positive vorticity over the warmer 
land surface and merged together (see Supplementary Figs. S8-S9 
online). Analysis from the experimental results also indicates that, 
the increase of LST (sensitivity experiments L3 and R3) leads to a 
convergence of moisture flux and latent heat at the low- middle atmo- 
sphere and an increase of PBL height over the warmer land surface, 
which favors the development of deep convection and thus supplies 
enough moisture and energy for the super typhoon (see 
Supplementary Figs. S10-S12 online). We can see more clearly the 
thermal-dynamical interactions between a strong typhoon and the 



atmosphere over the warmer land surface by examining the tem- 
poral-spatial variations of the four critical meteorological parameters 
for several strong typhoons: the apparent heat source (Q x ), apparent 
moisture sink (Q 2 ), and the vertical motion (see Supplementary Figs. 
S13-S18 online). The high value centers of Qi, Q 2 and vertical velo- 
city of the super typhoons appear to move toward those over warmer 
land surface. We found their interactions can be best sketched with a 
conceptual model illustrated by Fig. 6. During the typhoon season, 
the southwesterly EASM constantly transports water vapor from the 
southwestern sea to the warm land surface (indicated by the small 
white arrows) and the typhoon center (indicated by the large white 
arrows). There may be a convergence of moisture over the land where 
the moisture transported by EASM meets that transported by the 
detrained flows of the super typhoon (indicated by large yellow 
arrows). Meanwhile, warm land surface heats up the lower atmo- 
sphere, inducing stronger atmospheric convections. With the rise of 



SCIENTIFIC REPORTS | 3 : 1522 | DOI: 1 0.1 038/srep0 1522 



6 




Figure 6 | Cartoon of a conceptual model of the thermodynamic interactions between warmer land surface and a super typhoon approaching the coasts 
of China. The white and yellow arrows denote the moisture transportation by East Asian Summer Monsoon (EASM) and the detrained flows of the super 
typhoon, respectively. 



warm and moist surface air into the middle or upper troposphere, a 
large amount of latent heat is released, forming a warm core with 
anomalous Ch and Q 2 . To balance the upward motions, downward 
motions are also forced in the vicinity of the warm core. The intense 
heat and moisture over the warm land surface can feed typhoons in 
the same way as over the warm ocean surface, providing great incen- 
tives for typhoons to move toward the warmer land (see 
Supplementary Figs. S8-S18 online). 

Discussion 

It is worthy to note that the impact of warmer land surface on the 
super typhoon tracks is largely depended on the available moisture of 
the lower atmosphere that is necessary for the development of deep 
convection. In the southeastern China, the EASM brings plenty of 
water vapor from the ocean, which provides the necessary moisture 
for deep convection over the warmer land surface. In the regions 
where there is no moisture available, the impacts of warmer land 
surface on the super typhoon tracks may be small or negligible. 
Therefore, the applicability or accuracy of the LST predictor for 
typhoon tracks may be regionally dependent and a thorough invest- 
igation should be taken before using the LST as a TC track predictor 
for any other region. 

Since land and air temperatures are intimately coupled, we also 
examined if air temperature can serve as a reliable indicator to 
typhoon landfall tracks. Among all cases investigated, we found no 
significant correlation between landfall locations and air temperature 
over these locations or their adjacent areas. Higher air temperature 
over land surface is often linked to lower relative humidity and less 
large-scale ascent or deep convection. The dry circulation anomalies 
inhibit tropical cyclones from approaching 26 . Further analysis indi- 
cates that there is a good relationship between landfall locations of 
super typhoons and soil-air temperature differences (see Supple- 
mentary Figs. S19-S21 online). The reason could be that, the large 
soil -air temperature differences, along with the abundant moisture 
supplied by the EASM, are conducive to the large-scale ascent and 
deep convection required for the development and maintenance of 
tropical cyclones. 

It was estimated that half of the world's population will be living 
within 100 km of the coast by 2030 27 , which makes accurate predic- 
tion of the landfalls of super typhoons and hurricanes ever more 
important when making life-or-death decisions. It has been recog- 
nized that the mechanism to determine the movement of tropical 



cyclones is very complicated and a large part of it remains 
unresolved, especially for super typhoons. Although further invest- 
igation is still needed, our finding that super typhoons tend to make 
landfalls toward warmer land could be practically valuable to better 
cope with such potential natural disasters. Complementing the cur- 
rent steering flow theory, we propose the use of land surface temper- 
ature as a new indicator for improving landfall forecasts of super 
typhoons over southeastern China. Our results, while focusing on 
China, offer a new perspective on understanding the dynamic beha- 
viors of the most powerful windstorms on Earth. 

Methods 

The land surface temperature used in this study is obtained from the CMA archive of 
daily- mean land surface temperature measured at 753 stations over China (which has 
denser distribution in coastal regions of China) and interpolated to regular grid points 
with a resolution of 0.5°x 0.5° using the Cressman method 28 . A 9-point smoothing 
was applied 30 times to get the composite land surface temperature field at the time of 
2 days before landfalls of each super typhoon in Fig. 1, and three times to get the 
individual land surface temperature field in Fig. 2. The result for each grid point after 
the 9-point smoothing is a weighted average of the grid point plus the 8 surrounding 
points with the weights of 1 .0 (center point), 0.5 (the points at each side and above and 
below), and 0.3 (corner points). Any missing data points are not included in the sum; 
points beyond the grid boundary are considered to be missing. The location of highest 
land surface temperature 1-4 days before the landfall of a typhoon is identified by a 
searching algorithm within a coastal area of 30°X30°centered at typhoon center. 

The reanalysis data from the National Centers for Environmental Prediction 
(NCEP) of the United States National Oceanic and Atmospheric Administration (US 
NOAA), available 4 times a day at 0000 UTC, 0600 UTC, 1200 UTC, and 1800 UTC 
with horizontal resolution of 2.5° X 2.5° (before 2000) or 1°X 1° (on/after 2000), are 
used to derive the daily- mean 500-hPa potential height, the steering flow, the vertical 
velocity, the apparent heat source (Q x ), the apparent moisture sink (Q 2 ), and the 
horizontal moisture flux (F). The steering flow is obtained following the traditional 
method 8 ' 29 , which is a height- and area-averaged large-scale flow between 925 hPa 
and 300 hPa and over a ring- area with 3—8 grids away from the typhoon center. Q x , 
Q 2 and F are calculated using the following formulas: 

Q 2 = -I§ = + V#+ "f] =Q 21 + Q 22 + Q 23 . (2) 

F=^Vq. (3) 

See the Supplementary Methods for detailed descriptions of each variable in the above 
formulas. 

To see the sensitivity of the landfalling track of Longwang (2005) to the land surface 
temperature, we design a set of experiments by imposing different temperature 
increment At on the right or left hand side of the observed landing location, and 
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compare the predicted tracks with that by a control run without any change of land 
surface temperature. The following is a list of description for each experiment: 

Rl: At = 5 K, on the right; 
R2: At = 10 K, on the right; 
R3: At = 15 K, on the right; 
LI: At = 5 K, on the left; 
L2: At = 10 K, on the left; 
L3: At = 15 K, on the left; 

CNTR: control run, i.e., no change of land surface temperature. 

Fig. S4 displays the distribution of land surface temperature after an increment 
At = 5 K is imposed on the right (Rl) or left (LI) hand side of the observed landfalling 
location. 

Two nested domains with horizontal resolutions of 60 km and 20 km, respectively, 
are employed in all the experiments, with the central location of (22.2° N, 124.0° E). 
There 36 layers in the vertical with the top at 10 hPa. Two sets of experiments, 
initializing at 0000 UTC 30 Sept. (about 3 days before landing) and at 0000 UTC 1 
Oct. 2005 (about 2 days before landing) respectively, are performed. The initial 
conditions (IC) and boundary conditions (BC) for the first guess are generated from 
the NCEP Final Analysis dataset of the Global Forecast System with 1 ° X 1 ° horizontal 
resolution and 6-h interval. See Supplementary Methods for detailed descriptions of 
the physical schemes employed in the experiments. 
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